arXivil504.02402vl [hep-ph] 9 Apr 2015 


1 


SISSA 18/2015/FISI 
IPMU15-0042 


PREDICTIONS FOR THE DIRAC CP VIOLATION PHASE IN THE 
NEUTRINO MIXING MATRIX 

S. T. petcovQ 

SISSA/INFN, Trieste, Italy, and 
Kavli IPMU (WPI), University of Tokyo, Tokyo, Japan 

I. GIRARDI and A. V. TITOV 
SISSA/INFN, Trieste, Italy 


Using the fact that the neutrino mixing matrix U = utu,/, where Ue and result from 
the diagonalisation of the charged lepton and neutrino mass matrices, we analyse the 
predictions based on the sum rules which the Dirac phase 6 present in U satisfies when 
Uu has a form dictated by, or associated with, discrete flavour symmetries and Ue has a 
“minimal” form (in terms of angles and phases it contains) that can provide the requisite 
corrections to Ui>, so that the reactor, atmospheric and solar neutrino mixing angles ^13, 

623 and 9i2 have values compatible with the current data. 
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1 . Introduction 

One of the major goals of the future experimental studies in neutrino physics is the 
searches for CP violation (CPV) effects in neutrino oscillations (see, e.g., Refs. 1, 

2, 3, 4). It is part of a more general and ambitious program of research aiming to 
determine the status of the CP symmetry in the lepton sector. 

In the case of the reference 3-neutrino mixing scheme, CPV effects in the flavour 
neutrino oscillations, i.e., a difference between the probabilities of z/; -A vy and 
S'/ vy oscillations in vacuum!^, P{i'i -A vy) and P(i>/ -A uy), IV = e,^,T, can 
be caused, as is well known, by the Dirac phase present in the Pontecorvo, Maki, 
Nakagawa and Sakata (PMNS) neutrino mixing matrix U. If the neutrinos with 
dehnite masses Vi, i = 1,2,3, are Majorana particles, the 3-neutrino mixing matrix 
contains two additional Majorana CPV phased. However, the flavour neutrino 
oscillation probabilities P(vi -A vy) and P{Di -A vy), 1,1' = e,ji,T, do not depend 
on the Majorana phase^^. Our interest in the CPV phases present in the neutrino 
mixing matrix is stimulated also by the intriguing possibility that the Dirac phase 
and/or the Majorana phases in U can provide the CP violation necessary for the 
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generation of the observed baryon asymmetry of the Universel^ (see also, e.g., 
Refs. 10, 11). 

In the framework of the reference 3-flavour neutrino mixing we will consider, 
the PMNS neutrino mixing matrix is always given hy U = UlU^, where Ue and 
Uu are 3x3 unitary matrices originating from the diagonalisation of the charged 
lepton and the neutrino (Majorana) mass terms. We will suppose in what follows 
that Ui, has a form which is dictated by, or associated with, symmetries (see, e.g.. 
Refs. 12, 13). In the present article we consider the following symmetry forms of 
Uv'- i) tri-bimaximal (TBM)li2®^, ii) bimaximal (BM) (or corresponding to the 
conservation of the lepton charge L' = L^, — (LC))^^^^^, iii) golden ratio 

type A (GRA)l 2 il 22 |^ golden ratio type B (GRB)I^, and v) hexagonal (HG)I^. 
For all these symmetry forms U^, can be written as 

U, = ^-1 U, Qo = R23{0^3)Rl2i&12)Qo , ( 1 ) 

where i? 23 ( 023 ) -^ 12 (^ 12 ) orthogonal matrices describing rotations in the 

2-3 and 1-2 planes, respectively, and and Qq are diagonal phase matrices each 
containing two phases. The phases in the matrix Qq give contribution to the Ma¬ 
jorana phases in the PMNS matrix. The symmetry forms of C/^ of interest, TBM, 
BM (LG), GRA, GRB and HG, are characterised by the same values of the angles 
9iq = 0 and 6*23 = but correspond to different fixed values of the angle 05(2 

thus of siv? 9^2, namely, to i) sin^ 0('2 = 1/3 (TBM), ii) sin^ 6 IJ '2 = 1/2 (BM (LC)), 
iii) sin^ 9^2 = (2-l-r)“^ = 0.276 (GRA), r being the golden ratio, r = {l-\-y/h)/2, iv) 
sm^ 9^2 = (3 —r)/4 = 0.345 (GRB), and v) sin^ 0/2 = 1/4 (HG). The best fit values 
(b.f.v.) and Icr errors of the three corresponding neutrino mixing parameters in the 
standard parametrisation of the PMNS matrix®, which we will employ, read^^ 

sin^ 012 = 0.308;°;°!^ , ( 2 ) 

sin^ 013 = 0.0234^^^9 > (3) 

sin2 023 = 0.437tO:°i , (4) 

where the quoted values correspond to neutrino mass spectrum with normal or¬ 
dering (NO); the values for spectrum with inverted ordering (10) found inl^ differ 
insignificantly. The minimal form of Ue of interest that can provide the requisite 
corrections to fA, so that the neutrino mixing angles 0 i 3 , 023 and 0 i 2 have values 
compatible with the current data, including a possible sizeable deviation of 023 from 
7 r/ 4 , includes a product of two orthogonal matrices describing rotations in the 2-3 
and 1-2 planes^^, i? 23 ( 6 * 23 ) 1 ^ 12 (^ 12 ): ^23 ^12 being two (real) anglesj^ This 

leads to the following parametrisation of the PMNS matrix U: 

U = Ri2 ( 0 J 2 ) R23 {9^23) R23 (9^23) Ri2 { 9 ’( 2 ) Qo , ( 5 ) 


For a detailed discussion of alternative possibilities see Ref. 27. 
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where = diag (l, e e , and 6*23 = — ’’"/d. Equation M can be recast in the 
fornJ^ 


U = Rl2{0t2M^)R23{d23) Rl2(e'^2) Q , ( 6 ) 


where we have defined $ = diag (l,e®'^, l), (j) being a CPV phase, 623 is a function 
of 0 I 3 , sin^ 023 = 1/2 — sin 0|3 cos 0|3 cos(a; — -ijj), and Q is a diagonal phase matrix. 
The phases in Q give contributions to the Majorana phases in the PMNS matrix. 
The angle 023, however, can be expressed in terms of the angles 023 and 0 i 3 of the 
PMNS matrix: 


sin^ 023 = 


\U>.3\" 

l-\Ue3f 


sin^ 023 - sin^ 0i3 
1 — sin^ 013 


(7) 


and the value of 823 is fixed by the values of 623 and 0 i 3 . 


2. Predicting the Dirac Phase in the PMNS Matrix 

In the scheme under discussion the four observables 0 i 2 , 023, 0 i 3 and the Dirac 
phase S in the PMNS matrix are functions of three parameters ^ 12 , ^23 and (p. As a 
consequence, the Dirac phase <5 can be expressed as a function of the three PMNS 
angles 0 i 2 , 023 and 0 i 3 , leading to a new “sum rule” relating S and 0 i 2 , 623 and 0 i 3 . 
Within the approach employed this sum rule is exact. Its explicit form depends 
on the symmetry form of the matrix [/j^, i.e., on the value of the angle 0 / 2 - For 
arbitrary fixed value of 0/2 the sum rule of interest read^^ 

cos5 = — ^^n^/3 _ _ \oq^29'P2 + (sin^ 0 i 2 — cos^ 0 / 2 ) (l — cot^ 023 sin^ ^ 13 )] • 
Sin 26 'i 2 sin c'13 

( 8 ) 

A similar sum rule can be derived for the phase (p 

In Refs. 28, 30 we have derived predictions for cos <5, <5 and the rephasing invariant 
Jcp, which controls the magnitude of the CPV effects in neutrino oscillation^^, 
using the sum rule in eq. Q and the measured values of the lepton mixing angles 
^ 12 , ^13 and 023. In the present article we first summarise the predictions for these 
observables obtained in Refs. 28, 30 in a simplified analysis employing the best fit 
values (b.f.v.) and the 3 ct allowed ranges of the three relevant neutrino mixing 
parameters sin^ 0 i 2 , sin^ 0 i 3 and sin^ 023. This is followed by a summary of the 
results of the statistical analysis of the predictions performed in Ref. 30, which is 
based on i) the current, and most importantly, ii) the prospective, uncertainties in 
the measured values of sin^ 0 i 2 , sin^ 0 i 3 and sin^ 023 . 

We note first that the predicted values of cos <5 vary significantly with the sym¬ 
metry form of Ui,^. For the best fit values of sin^ 0 i 2 = 0.308, sin^ 0 i 3 = 0.0234 
and sin^ 023 = 0 . 437 found inl^, for instance, we get cos<5 = (—0.0906), (—1.16), 
0.275, (-0.169) and 0.445 for the TBM, BM (LC), GRA, GRB and HG forms, 
respectively. For the TBM, GRA, GRB and HG forms these values correspond 
to (5 = ±95.2°, ±74.0°, ±99.7°, ±63.6°, respectively. The unphysical value of cos <5 
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in the BM (LC) case is a reflection of the fact that the scheme under discussion 
with BM (LC) form of the matrix ll^, does not provide a good description of the 
current data on 0 i 2 , 023 and 013 Physical values of cosi5 can be obtained, for 
instance, for the b.f.v. of sin^ 0 i 3 and sin^ 023 if sin^ 012 has a larger value!^ for, 
e.g., sin^0i2 = 0.34 allowed at 2a by the current data, we have cos 5 = —0.943, 
corresponding to 5 = ±160.6°. Similarly, for sin^ 0i2 = 0.32, sin^ 023 = 0.41 and 
sin 013 = 0.158 we haveP^ cosi5 = —0.978, <5 = ±168.1°. 

The results quoted above implj^ that the measurement of cos i5 can allow to 
distinguish between the different symmetry forms of Uv, provided 0 i 2 , 0 i 3 and 023 
are known with a sufficiently good precision. Even determining the sign of cos 5 will 
be sufficient to eliminate some of the possible symmetry forms of Un. 

It was also found in^^ that the sum rule predictions for cos 6 exhibit strong 
dependence on the value of sin^ 0 i 2 when the latter is varied in its ia experimentally 
allowed range^^ (0.259-0.359). The predictions for cos 5 change significantly not 
only in magnitude, but also the sign of cos 5 changes in the TBM, GRA and GRB 
case^^. In the case of 0^3 = 0, for instance, we get for the TBM form of Ui, for 
the three values of sin^ 0i2 = 0.308, 0.259 and 0.359: cosJ = (—0.114), (—0.469) 
and 0.221, thus cos (5 = 0 is allowed for a certain value of sin^0i2. For the GRA 
and GRB forms of IJ,y we have, respectively, cos(5 = 0.289, (—0.044), 0.609, and 
cos(5 = (—0.200), —0.559, 0.138. Similarly, for the HG form we find for the three 
values of sin^0i2: cos(5 = 0.476, 0.153, 0.789. 

In what concerns the dependence of the sum rule predictions for cos S when 
sin^ 023 is varied in its 3tT allowed interval, 0.374 < sin^ 023 < 0.626, the results we 
obtained for sin^ 023 = 0.374 and sin^ 023 = 0.437, setting sin^ 0i2 and sin^ 0i3 to 
their best fit values, do not differ significantly. However, the differences between 
the predictions for cos <5 obtained for sin^ 023 = 0.437 and for sin^ 023 = 0.626 are 
rather largely — they differ by the factors of 2.05, 1.25, 1.77 and 1.32 in the TBM, 
GRA, GRB and HG cases, respectively. 

Similar analysis can be performed for the predictions for the cosine of the phase 
(j) which in many theoretical models serves as a “source” for the Dirac phase S. 
The phase (j) is related to, but does not coincide with, the Dirac phase S. This leads 
to the confusing identification of (j) with 6 : the sum rules satisfied by cos (j) and cos <5 
differ significantly!^. Correspondingly, the predicted values of cos (j) and cos S in the 
cases of the TBM, GRA, GRB and HG symmetry forms of Ui, considered by us also 
differ significantly. This conclusion is not valid for the BM (LC) form: for this form 
the sum rules predictions for cos (p and cos S are rather similar!^. 

We next present results of the statistical analysis of the predictions for S, cos <5 
and the rephasing invariant Jqp performed in Ref. 30 in the cases of the TBM, 
BM (LC), GRA, GRB and HG symmetry forms of the matrix U^. In this analysis 
the latest results on sin^0i2, sin^0i3, sin^ 023 and S, obtained in the global anal¬ 
ysis of the neutrino oscillation data performed in^^ were used as input. The aim 
was to derive the allowed ranges for cos 5 and Jcp, predicted on the basis of the 
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current data on the neutrino mixing parameters for each of the symmetry forms of 
Un considered. For this purpose the function was constructed in the following 
waj® X^i{xi}) = with Xi = {sin^ 6 »i 2 , sin^ 6 »i 3 , sin^ 6»23, i5}- The func¬ 

tions Xi have been extracted from the 1-dimensional projections given in Ref. 25 
and, thus the correlations between the oscillation parameters have been neglected. 
This approximation is sufhciently precise since it allows to reproduce the contours in 
the planes (sin^ 023 , ^), (sin^ 0 i 3 , 5) and (sin^ 023 , sin^ 0 i 3 ), given in^^, with a rather 
high accuracy. We calculated x^(cos(5) by marginalising over sin^ 0 i 3 and 

sin^ 023 for a fixed value of cos 6 . Given the global fit results, the likelihood function, 

, (9) 

represents the most probable value of cos (5 for each of the considered symmetry 
forms of Ui,. The na confidence level (C.L.) region corresponds to the interval of 
values of cos<5 in which L{cosd) > L{x^ = Xmin) ' = n^), where Xmin is ^he 

value of x^ in the minimum. 

In Fig. we show the likelihood function versus cos 5 for NO neutrino mass 
spectrum from Ref. 30. The results shown are obtained by marginalising over 
all the other relevant parameters of the scheme considered. The dependence of 
the likelihood function on cos 6 in the case of 10 neutrino mass spectrum differs 
little from that shown in Fig. As can be observed in Fig. a rather precise 
measurement of cos S would allow to distinguish between the different symmetry 
forms of Ui, considered by us. For the TBM and GRB forms there is a significant 
overlap of the corresponding likelihood functions. The same observation is valid 
for the GRA and HG forms. However, the overlap of the likelihood functions of 
these two groups of symmetry forms occurs only at 3cr level in a very small interval 
of values of cosi5. This implies that in order to distinguish between TBM/GRB, 
GRA/HG and BM (LC) symmetry forms, a not very demanding measurement (in 
terms of accuracy) of cos (5 might be sufficient. The value of the non-normalised 
likelihood function at the maximum in Fig. is equal to exp(—x^j^/S), which 
allows us to make conclusions about the compatibility of the symmetry schemes 
considered with the current global data. The results of this analysis for cos 5 are 
summarised in Table [T] 

We have also performed in Ref. 30 a similar statistical analysis of the predictions 
for the rephasing invariant Jcp in the cases of the TBM, BM (LC), GRA, GRB 
and HG symmetry forms of the matrix C/jx considered. In this analysis we used as 
input the latest results on sin^ 0 i 2 , sin^ 0 i 3 , sin^ 023 and S, obtained in the global 
analysis of the neutrino oscillation data performed in Ref. 25, and minimised x^ 
for a fixed value of Jcp. The obtained b.f.v. and 3a ranges are given in Table 
We have found, in particular, that the CP-conserving value of Jcp = 0 is excluded 
in the cases of the TBM, GRA, GRB and HG neutrino mixing symmetry forms, 
respectively, at approximately 5a, 4a, 4a and 3 ct C.L. with respect to the C.L. of the 
corresponding best fit value. These results reflect the predictions we have obtained 
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Fig. 1. The likelihood function versus cos <5 for NO neutrino mass spectrum after marginalising 
over sin^ @13 and sin^ 023 1 for the TBM, BM (LC), GRA, GRB and HG symmetry forms of the 
mixing matrix Uiy (see text for further details). (From Ref. 30.) 


for S, more specifically, the C.L. at which the CP-conserving values of <5 = 0 
TT, are excluded in the discussed cases. We found also that the 3a allowed intervals 
of values of <5 and Jcp are rather narrow for all the symmetry forms considered, 
except for the BM (LC) form. More specifically, for the TBM, GRA, GRB and 
HG symmetry forms we have obtained at Scr: 0.020 < |Jcp| < 0.039. For the 
b.f.v. of Jcp we have found, respectively: Jcp = (—0.034), (—0.033), (—0.034), 
and (—0.031). Our results indicate that distinguishing between the TBM, GRA, 
GRB and HG symmetry forms of the neutrino mixing would require extremely high 
precision measurement of the Jcp factor^. 

In Fig. we present the likelihood function versus cos <5 within the 
Gaussian approximation, i.e., using Xq = with yi = 

{sin^ 012 , sin^ 6 * 13 , sin^ 023 }, where we used the current b.f.v. (j/j) of the mixing 
angles for NO neutrino mass spectrum given in Ref. 25 and the prospective la 
uncertainties {ay.) in the determination of sin^ 0 i 2 (0.7% from JUNOl^, sin^ 0 i 3 
(3% derived from an expected error on sin^ 20 i 3 of 3% from Daya Bay, see Refs. 4, 
32) and sin^ 023 (5% derived from the potential sensitivity of NOvA and T2K on 
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Table 1. Best fit values (b.f.v.) of Jqp and cos <5 and corresponding 3 (j ranges (found fixing 
— Xmin “ setup using the data froml^ for NO neutrino mass spectrum. (Prom 

Ref. 30, where results for 10 spectrum are also given.) 


Scheme 

Jcp/10 2 (b.f.v.) 

•Icp/lO ^ (3cr range) 

COS 6 (b.f.v.) 

COS 5 (3(J range) 

TBM 

-3.4 

[-3.8,-2.8] U [3.1, 3.6] 

-0.07 

[- 0 . 47 , 0 . 21 ] 

BM (LC) 

-0.5 

[-2.6, 2.1] 

-0.99 

[-1.00,-0.72] 

GRA 

-3.3 

[-3.7,-2.7] U [3.0, 3.5] 

0.25 

[-0.08, 0.69] 

GRB 

-3.4 

[-3.9,-2.6] U [3.1, 3.6] 

-0.15 

[-0.57, 0.13] 

HG 

-3.1 

[-3.5,-2.0] U [2.6, 3.4] 

0.47 

[ 0.16, 0.80] 


sin^ 2023 of 2%, see Ref. 4, this sensitivity can be also achieved in future neutrino 
facilities as T2HKI^. The BM (LC) case is very sensitive to the b.f.v. of sin^ 0i2 
and sin^ 023 and is disfavoured at more than 2 a for the current b.f.v. found in^^. 
This case might turn out to be compatible with the data for larger (smaller) mea¬ 
sured values of sin^ 0i2 (sin^ S 23 ), as can be seen from Fig. which was obtained 
for sin^ 0i2 = 0.332 (the best fit values of the two other mixing angles being kept 
intact). With the increase of the value of sin^ 023 the BM (LC) form becomes in¬ 
creasingly disfavoured, while the TBM/GRB (GRA/HG) predictions for cos <5 are 
shifted somewhat to the left (right) with respect to those shown in Fig. For, e.g., 
the best fit values of sin^ 0i2 = 0.304, sin^ 0i3 = 0.0219 and sin^ 623 = 0.579, found 
in Ref. 34 for 10 neutrino mass spectrum, these shifts in cos 6 are approximately by 
0 . 1 . 

The measurement of sin^ 0i2, sin^ 0i3 and sin^ 023 with the quoted precision will 
open up the possibility to distinguish between the BM (LC), TBM/GRB, GRA 
and HG forms of U,^. Distinguishing between the TBM and GRB forms would 
require relatively high precision measurement of cos S. Assuming that | cos (5| < 0.93, 
which means for 76% of values of S, the error on S, AS, for an error on cos 5, 
A(cos^) = 0.10(0.08), does not exceed AJ < A(cos^)/-\/l — 0.93^ = 16° (12°). 
This accuracy is planned to be reached in the future neutrino experiments like 
T2HK (ESSi^SB)®. Therefore a measurement of cos (5 in the quoted range will 
allow one to distinguish between the TBM/GRB, BM (LG) and GRA/HG forms at 
approximately 3cr G.L., if the precision achieved on sin^ 0i2, sin^ 0i3 and sin^ 023 is 
the same as in Figs. and 

3. Summary and Conclusions 

In conclusions, we have derived inl^Sl the ranges of the predicted values of cos S and 
Jcp for the TBM, BM (LG), GRA, GRB and HG symmetry forms of from a 
statistical analysis using the sum rule in eq. (|^ obtained in^^ and the current global 
neutrino oscillation data^^. The results of this analysis are summarised in Table 
and in Fig.[^ We found, in particular, that in the TBM, GRA, GRB and HG cases, 
the best fit values of Jcp lie in the narrow interval (—0.034) < Jcp < (—0.031), 
while at 3cr we have 0.020 < |Jcp| < 0.039. The predictions for 5, cos 5 and Jcp 







Fig. 2. The same as in Fig. but using the prospective Icr uncertainties in the determination 
of sin^ 0 i 2 , sin^ 013 and sin^ 623 within the Gaussian approximation. The three neutrino mixing 
parameters are fixed to their current best fit values (i.e., sin^ 0 i 2 = 0.308, etc.). See text for 
further details. (From Ref. 30.) 


in the case of the BM (LC) symmetry form of U,^, as the results of the statistical 
analysis performed by us showed, differ significantly: the best fit value of 5 = tt, 
and, correspondingly, of Jqp — 0. For the Scr range in the case of NO (10) neutrino 
mass spectrum we find: —0.026 (—0.025) < Jcp < 0.021 (0.023), i.e., it includes a 
sub-interval of values centred on zero, which does not overlap with the 3 (t allowed 
intervals of values of Jcp, corresponding to the TBM, GRA, GRB and HG symmetry 
forms of U,^. 

Finally, we have derived in® predictions for cos <5 using the prospective Ict uncer¬ 
tainties in the determination of sin^ 0 i 2 , sin^ ^13 and sin^ 023 respectively in JUNO, 
Daya Bay and accelerator and atmospheric neutrino experiments (Figs. and |^. 
The results thus obtained show that i) the measurement of the sign of cos <5 will 
allow to distinguish between the TBM/GRB, BM (LC) and GRA/HG forms of U^, 
ii) for a best fit value of cos <5 = —1 (—0.1) distinguishing at 3cr between the BM 
(TBM/GRB) and the other forms of Ui, would be possible if cos <5 is measured with 
Ict uncertainty of 0.3 (0.1). 
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Fig. 3. The same as in Fig.[^ but using sin^ 9x2 = 0.332. (From Ref. 30.) 

The results obtained in the studies performed in Refs. 28, 30 show, in particular, 
that the experimental measurement of the Dirac phase 5 of the PMNS neutrino 
mixing matrix in the future neutrino experiments, combined with the data on the 
neutrino mixing angles can provide unique information about the possible discrete 
symmetry origin of the observed pattern of neutrino mixing. 
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